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I The automation of the traditional Painleve test in Mathematica is discussed. The 

package PainleveTest .m allows for the testing of polynomial systems of nonlinear 
ordinary and partial differential equations which may be parameterized by arbitrary 
^0 I functions (or constants). Except where limited by memory, there is no restriction 

pH . on the number of independent or dependent variables. The package is quite robust 

in determining all the possible dominant behaviors of the Laurent series solutions 
of the differential equation. The omission of valid dominant behaviors is a common 
problem in many implementations of the Painleve test, and these omissions often 
. lead to erroneous results. Finally, our package is compared with the other available 

^ ' implementations of the Painleve test. 

o 

O ; 1 Introduction 

, Completely integrable nonlinear partial differential equations (PDEs) have remarkable 

I properties, such as infinitely many generalized symmetries, infinitely many conservation 

laws, the Painleve property, Backlund and Darboux transformations, bilinear forms, and 
Lax pairs (cf. [2,11,24,25]). Completely integrable equations model physically interesting 
wave phenomena in reaction-diffusion systems, population and molecular dynamics, non- 

■ linear networks, chemical reactions, and waves in material science. By investigating the 
^ . complete integrability of a nonlinear PDE, one gains important insight into the structure 

^ I of the equation and the nature of its solutions. 

Broadly speaking, Painleve analysis is the study of the singularity structure of differen- 
tial equations. Specifically, a differential equation is said to have the Painleve property if all 
the movable singularities of all its solutions are poles. There is strong evidence [48,50,51] 
that integrability is closely related to the singularity structure of the solutions of a differ- 
ential equation (cf. [33,38]). For instance, dense branching of solutions around movable 
singularities has been shown to indicate nonintegrability [49]. 

At the turn of the nineteenth-century, Painleve [30] and his colleagues classified all the 
rational second-order ODEs for which all the solutions are single- valued around all movable 
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singularities. Equations possessing this property could either be solved in terms of known 
functions or transformed into one of the six Painleve equations whose solutions define the 
Painleve transcendents. The Painleve transcendents cannot be expressed in terms of the 
classical transcendental functions, except for special values of their parameters [19]. 

The complex singularity structure of solutions was first used by Kovalevskaya in 1889 
to identify a new integrable system of equations for the motion for a rotating top (cf. [14, 
38]). Ninety years later, Ablowitz, Ramani and Segur (ARS) [2, 3] and McLeod and 
Olver [27] formulated the Painleve conjecture which gives a useful necessary condition 
for determining whether a PDE is solvable using the Inverse Scattering Transform (1ST) 
method. Specifically, the Painleve conjecture asserts that every nonlinear ODE obtained 
by an exact reduction of a nonlinear PDE solvable by the IST-method has the Painleve 
property. While necessary, the condition is not sufficient; in general, most PDEs do 
not have exact reductions to nonlinear ODEs and therefore satisfy the conjecture by 
default [41]. Weiss, Tabor and Carnevale (WTC) [44] proposed an algorithm for testing 
PDEs directly (which is analogous to the ARS algorithm for testing ODEs) . For a thorough 
discussion of the traditional Painleve property, see [1,8,10,13,18,28,31,33,38,39]. 

There are numerous methods for solving completely integrable nonlinear PDEs, for 
instance by explicit transformations into linear equations or by using the IST-method [11]. 
Recently, progress has been made using Mathematica and Maple in applying the IST- 
method to difficult equations, including the Camassa-Holm equation [21]. While there 
is as yet no systematic way to determine if a differential equation is solvable using the 
IST-method [27], having the Painleve property is a strong indicator that it will be. 

There are several implementations of the Painleve test in various computer algebra 
systems, including Reduce, Macsyma, Maple and Mathematica. The implementations 
described in [34, 35, 37] are limited to ODEs, while the implementations discussed in [16, 
45-47] allow the testing of PDEs directly using the WTC algorithm. The implementation 
for PDEs written in Mathematica by Hereman et al. [16] is limited to two independent 
variables {x and t) and is unable to find all the dominant behaviors in systems with 
undetermined exponents ai (as is the case with the Hirota-Satsuma system). Our package 
PainleveTest .m [4] written in Mathematica syntax, allows the testing of polynomial PDEs 
(and ODEs) with no limitation on the number of differential equations or the number of 
independent variables (except where limited by memory). Our implementation also allows 
the testing of differential equations that have undetermined dominant exponents Oi and 
that are parameterized by arbitrary functions (or constants). The implementations for 
PDEs written in Maple by Xu and Li [45-47] were written after the one presented in this 
paper and are comparable to our implementation. 

The paper is organized as follows: in Section[2lwe review the basics of Painleve analysis. 
Section |51 discusses the WTC algorithm for testing PDEs and uses the Korteweg-de Vries 
(KdV) equation and the Hirota-Satsuma system of coupled KdV (cKdV) equations to 
show the subtleties of the algorithm. We detail the algorithms to determine the dominant 
behavior, resonances, and constants of integration using a generalized system of coupled 
nonlinear Schrodinger (NLS) equations in Section |11 Additional examples are presented 
in Section [3 to illustrate the capabilities of the software. Section El compares our software 
package to other codes and briefly discusses the generalizations of the WTC algorithm. The 
use of the package PainleveTest .m [4] is shown in Section [7| We draw some conclusions 
and discuss the results in Section |H1 
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2 Painleve Analysis 

Consider a system of M polynomial differential equations, 

F,(u(z), u'(z), u"(z), . . . , u(™-) (z)) = 0, i = 1, 2, . . . , M, (2.1) 

where the dependent variable u(z) has components mi(z), . . . , um{z), the independent vari- 
able z has components zi, . . . , zn, and u('^)(z) = d^u{z)/{dzi^dz2^ ■ ■ ■ dz^) denotes the 
collection of mixed derivative terms of order k. Let m = Xli^i "^ii where rrii is the highest 
order in each equation. If there are any arbitrary coefficients (constants or analytic func- 
tions of z) parameterizing the system, we assume they are nonzero. For simplicity, in the 
examples we denote the components of u(z) by u(z), v{z),'w{z), . . . , and the components 
of z by X, y, . . . , t. 

A differential equation has the Painleve property if all the movable singularities of all its 
solutions are poles. A singularity is movable if it depends on the constants of integration 
of the differential equation. For instance, the Riccati equation, 

w'{z)+w'^{z) = 0, (2.2) 

has the general solution w{z) = l/{z — c), where c is the constant of integration. Hence, 
(|2.2|) has a movable simple pole at z = c because it depends on the constant of integration. 
Solutions of ODEs can have various kinds of singularities, including branch points and 
essential singularities; examples of the various types of singularities [23] are shown in 
Tabled As a general property, solutions of linear ODEs have only fixed singularities [19]. 

Simple fixed pole 

z'w' + 711 = ^ w{z)=c/z 

Simple movable pole 

w' + w'^ = => u){z) = l/{z - c) 

Movable algebraic branch point 

2ww' — 1 = =^ w{z) = \J z — c 

Movable logarithmic branch point 

w" + w'"^ = =^ w{z) = log(z — Cl) + C2 

Non-isolated movable essential singularity 

(1 + u?)w" + (1 - 2w)w''^ = w{z) = tan{ln(ci2; + C2)} 

Table 1. Examples of various types of singularities. 



In general, a function of several complex variables cannot have an isolated singular- 
ity [29]. For example, f{z) = 1/z has an isolated singularity at the point z = 0, but 
the function f{w, z) = 1/z of two complex variables, w = u + iv,z = x + iy, has a two- 
dimensional manifold of singularities, namely the points (n, v,0,0), in the four-dimensional 
space of these variables. Therefore, we will define a pole of a function of several complex 
variables as a point (ai, 02, ... , ajy), in whose neighborhood the function can be written 
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in the form /(z) = h(z)/g{z), where g and h are both analytic in a region containing 
(oi, . . . , Oat) in its interior, g{ai, . . . , a^) = 0, and /i(ai, . . . , oat) ^ 0. 

The WTC algorithm considers the singularity structure of the solutions around non- 
characteristic manifolds of the form g{z) = 0, where g{z) is an analytic function of z = 
{zi,Z2, ■ ■ ■ , zn) in a neighborhood of the manifold. Specifically, if the singularity manifold 
is determined by g{z) = and u(z) is a solution of the PDE, then one assumes a Laurent 
series solution 

oo 

n,(z) =g"'(z) i = l,2,...,M, (2.3) 

fc=0 

where the coefficients Ui^k{z) are analytic functions of z with Uifi{z) ^ in a neighborhood 
of the manifold and the ai are integers with at least one exponent < 0. The requirement 
that the manifold g{z) = is non-characteristic, ensures that the expansion (|2.3() is well 
defined in the sense of the Cauchy-Kovalevskaya theorem [41,43]. 

Substituting (|2.3j) into (|2.H) and equating coefficients of like powers of g{z) determines 
the possible values of and defines a recursion relation for Uj /^(z). The recursion relation 
is of the form 

QfeUfc = Gfc(uo,ui, . . . ,Ufe„i,g,z), = {ui^k,U2,k,- ■ ■ ,UM,k)^ 1 (2.4) 

where Qk is an M x M matrix and T denotes transpose. 

For (|2.H) to pass the Painleve test, the series p.3|) should have m — 1 arbitrary functions 
as required by the Cauchy-Kovalevskaya theorem (as g{z) is the m-th arbitrary function). 
If so, the Laurent series solution corresponds to the general solution of the equation [1]. 
The m — 1 arbitrary functions Ui^ki'^) occur when k is one of the roots of det{Qk)- These 
roots ri < r2 < • • • < are called resonances. The resonances are also equal to the Fuchs 
indices of the auxiliary equations of Darboux [7] . 

Since the WTC algorithm is unable to detect essential singularities, it is only a nec- 
essary condition for the PDE to have the Painleve property [6]. While rarely done in 
practice, sufficiency is proved by finding a transformation which linearizes the differen- 
tial equation, yields an auto-Backlund transformation, a Backlund transformation, or 
hodographic transformation [15] to another differential equation which has the Painleve 
property (see [8,23,33] for more information). 

3 Algorithm for the Painleve Test 

In this section, we outline the WTC algorithm for testing PDEs for the Painleve property. 
We discuss the Kruskal simplification and the Painleve test of ODEs after the three main 
steps are outlined. Each of these steps is illustrated using both the KdV equation and 
the cKdV equations due to Hirota and Satsuma. Details of the three main steps of the 
algorithm are postponed till Sectional 

Step 1 (Determine the dominant behavior). It is sufficient to substitute 

ni(z) =Xi<7°'(z), i = l,2,...,M, (3.1) 

where Xi is a constant, into 1)2. 1|) to determine the leading exponents Oj. In the resulting 
polynomial system, equating every two or more possible lowest exponents of g{z) in each 
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equation gives a linear system for a^. The linear system is then solved for and each 
solution branch is investigated. The traditional Painleve test requires that all the a, are 
integers and that at least one is negative. 

If any of the ai are non- integer in a given branch, then that branch of the algorithm 
terminates. A non-integer Oj implies that some solutions of (|2.H) have movable algebraic 
branch points. Often, a suitable change of variables in 1)2. ip can remove the algebraic 
branch point. An alternative approach is to use the "weak" Painleve test, which allows 
certain rational ai and resonances; see [13,18,32,33] for more information. 

If one or more ai remain undetermined, we assign integer values to the free so that 
every equation in (|2.1|) has at least two different terms with equal lowest exponents. 

For each solution ai, we substitute 

ni(z) =u,,o(z)5"'(z), i = l,2,...,M, (3.2) 

into (|2.H) . We then solve the (typically) nonlinear equation for Uifi{z), which is found 
by balancing the leading terms. By leading terms, we mean those terms with the lowest 
exponent of g{z). If any of the solutions contradict the assumption that Uifi{z) ^ 0, then 
that branch of the algorithm fails the Painleve test. 

If any of the Oj are non-integer, all the Oj are positive, or there is a contradiction with 
the assumption that Uj.o(z) ^ 0, then that branch of the algorithm terminates and does 
not pass the Painleve test for that branch. 

Step 2 (Determine the resonances). For each ai and Ui^o{z), W6 CcllculcltG th.G ^1 ^ * * * ^ f^m 
for which Ui^r{^) is an arbitrary function in (|2.3p . To do this, we substitute 

Ui{z) = u,fl{z)g'^^ (z) + ni,,(z)5"»+^(z) (3.3) 

into (|2.1|) . and keep only the lowest order terms in g{z) that are linear in Ui^r- This is done 
by computing the solutions for r of det(Qr) = 0, where the M x M matrix Qr satisfies 

QrUr = 0, Ur = {Ui^r U2,r ■■■ UM,r)'^ ■ (3.4) 

If any of the resonances are non-integer, then the Laurent series solutions of (|2.1|) have 
a movable algebraic branch point and the algorithm terminates. If r^ is not a positive 
integer, then the algorithm terminates; if ri = — l,r2 = • • • = = and m — 1 of 
the Uifi{z) found in Step 1 are arbitrary, then ()2.1j) passes the Painleve test. If ()2.1(l is 
parameterized, the values for ri < • • • < may depend on the parameters, and hence 
restrict the allowable values for the parameters. 

There is always a resonance r = — 1 which corresponds to the arbitrariness of g{z); as 
such, it is often called the universal resonance. When there are negative resonances other 
than r = —1, (or, more than one resonance equals —1) then the Laurent series solution 
is not the general solution and further analysis is needed to determine if (|2.H) passes the 
Painleve test. The perturbative Painleve approach, developed by Conte et al. [9], is one 
method for investigating negative resonances. 

Step 3 (Find the constants of integration and check compatibility conditions). For the 
system to possess the Painleve property, the arbitrariness of Ui,r{z) must be verified up to 
the highest resonance level. This is done by substituting 

Ui{z)=g»^iz)^Ui^k{z)g\z) (3.5) 

k=0 
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into (|2.H) . where is the largest positive integer resonance. 

For ()2.1() to have the Painleve property, the (M + 1) x M augmented matrix {Qk\Gk) 
must have rank M when k ^ r and rank M — s when k = r, where s is the algebraic 
multiplicity of r in det{Qr) = 0, 1 < k < r^, and Qk and are as defined in ()2.4() . 
If the augmented matrix {Qk\Gk) has the correct rank, solve the linear system p.4|) for 
ui^fc(z), . . . , UM,ki^) and use the results in the linear system at level k + 1. 

If the linear system (j2.4j) does not have a solution, then the Laurent series solution of 
1)2. 1(1 has a movable logarithmic branch point and the algorithm terminates. Often, when 
(|2.1|) is parameterized, carefully choosing the parameters will resolve the difference in the 
ranks of Qk and {Qk\Gk)- 

If the algorithm does not terminate, then the Laurent series solutions of 1)2.1(1 are free 
of movable algebraic or logarithmic branch points and ()2.1j) passes the Painleve test. 

The Painleve test of PDEs is quite cumbersome; in particular. Step 3 is lengthy and 
prone to error when done by hand. To simplify Step 3, Kruskal proposed a simplifica- 
tion which now bears his name. In the context of the WTC algorithm, it is sometimes 
called the Weiss-Kruskal simplification [20,23]. The manifold defined by g{z) = is non- 
characteristic, that means gzii'^) 7^ for some / on the manifold g{z) = 0. By the implicit 
function theorem, we can then locally solve g{z) = for zi, so that 

g{z) = zi- h{zi, ... , zi^i,zi^i, ...,zn), (3.6) 

for some arbitrary function h. Using ()3.6() greatly simplifies the computation of the con- 
stants of integration . . . , zi^i, . . . , zn)- However, with the Kruskal simplifica- 
tion one looses the ability to use the Weiss truncation method [42] to find a linearising 
transformation, an auto-Backlund transformation, or a Backlund transformation (see [8]). 
When testing ODEs, ()2.3j) must be replaced by 

oo 

u,{z)=g"'{z)Y,u^,k9H^), i = l,2,...,M, (3.7) 

k=0 

where the coefficients Ui^k are constants, g{z) = z — zq, and zq is an arbitrary constant. If 
z explicitly occurs in the ODE, then it is (automatically) replaced by g{z) + zq prior to 
Step 1 of the test. An example of the Painleve test of an ODE is given in Sectional 



3.1 The Korteweg-de Vries equation 

To illustrate the steps of the algorithm, let us examine the KdV equation [1], 

ut + 6uux + u^x = 0, (3.8) 

the most famous completely integrable PDE from soliton theory. Note that for simplicity, 
we use = Uxx---x = d^u/dx^ and gix = d^g/dx^ when i>2>. 
Substituting ()3.1|) into (|3.8|) gives 

ax{fftff""'+6x5x5'""'+5°"'[(«-l) ((«-2)5x+35fc) Sx+fiV]} = 0. (3.9) 

The lowest exponents of g{x,t) are a — 3 and 2a — 1. Equating these leading exponents 
gives a = —2. Substituting ()3.2|) . u{x,t) = uo{x,t)g~'^{x,t), into ()3.8() and requiring that 
the leading terms (in g~^{x,t)) balance, gives uo{x,t) = —2g'^{x,t). 
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Substituting (|3.3|) . u{x,t) = —2g1{x,t)g '^{x,t)+Ur{x,t)g'^ ^{x,t), into (|3.8|) and equat- 
ing the coefficients of the dominant terms (in g^~^{x,t)) that are hnear in Ur{x,t) gives 

(r - 6)(r - 4)(r + l)g^{x, tf = 0. (3.10) 

Assuming gx{x,t) ^ 0, the resonances of (|3.8() are ri = — l,r2 = 4 and = 6. 
We now substitute 

6 

u{x, t) = g~'^{x, t) ^ Uk{x, t)g^{x, t) 

= -2gl{x, t)g-\x, t) + ui{x, t)g-\x, *) + ••• + uq{x, t)g\x, t) (3.11) 

into (|3.8j) and group the terms of hke powers of g{x^ t). So, we will pull off the coefficients 
of g''^^{x,t) at level k. Equating the coefficients of g~^{x,t) to zero at level k = l, gives 
ui{x,t)gl{x,t) = 2gl{x,t)gxx{x,t). Setting tii(x, t) = 2fc(x, t), we get 



9t9l + 3ffxgL - 45'x5'3; 
at level k = 2. Similarly, at level k = 3, 



U2[x,t)= ^^^3 , (3.12) 



/ _ olgxt - gtQxQxx + '^gtx ~ ^9xgxx93x + glg^x rrno\ 

U3[X,l) — . ^O.iOj 

At level /c = r2 = 4, we find 

{Ui)t + Q{u3{uo)x + U2(u.l)x + Ui{u2)x + ?io(^i3)x} + ("uOsx = 0, (3.14) 

which is trivially satisfied upon substitution of the solutions oi uo{x,t), . . . ,us{x,t). There- 
fore, the compatibility condition at level A: = r2 = 4 is satisfied and U4^{x, t) is indeed arbi- 
trary. At level k = 5, U5{x,t) is unambiguously determined, but not shown due to length. 
Finally, the compatibility condition at level A; = = 6 is trivially satisfied when the 
solutions for UQ{x,t), . . . ,U3{x,t) and u^{x,t) are substituted into the recursion relation 
at that resonance level. 

Therefore, the Laurent series solution u{x, t) of (|3.(S)) in the neighborhood of g{x, t) = 
is free of algebraic and logarithmic movable branch points. Furthermore, since the Laurent 
series solution, 

oo 

u{x, t) = g-\x, t) Uk{x, t)g^{x, t), (3.15) 

fc=0 

has three arbitrary functions, g{x,t),U4{x,t), and UQ{x,t), (as required by the Cauchy- 
Kovalevskaya theorem since (|3.8|) is of third order) it is also the general solution. Hence, 
we conclude that (|3.8|) passes the Painleve test. 

The Weiss-Kruskal simplification uses g{x,t) = x — h{t). Consequently, gx = l,gxx = 
gsx = • • • = 0, and the Laurent series, 

oo 

u{x, t) = g-'^{x, t) u,^k{t) g'^ix, t), (3.16) 

fc=0 
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becomes 

1 



t) = --^—j-^ + ^h'{t) + Ui{t){x - h{t)f 



+ ^h"{t){x - h{t)f + U(^{t){x - h{t)f + • • • , (3.17) 
where h{t),U4{t) and UQ{t) are arbitrary. 

3.2 The Hirota-Satsuma system 

To show the subtleties in determining the dominant behavior, consider the cKdV equations 
due to Hirota and Satsuma [1] with real parameter a, 

ut = a{QuUx + U3x) - 2ufx, a > 0, 

(o.loj 

Vt = -3uVx - V3x. 

Again, we substitute ()3.1|) . u{x,t) = xig°'^{x,t) and v{x,t) = X2g°'^{x,t), into ()3.18() and 
pull off the lowest exponents of g{x,t). ^From the first equation, we get ai — 3,2ai — 1, 
and 2a2 — 1. ^From the second equation, we get 02 — 3 and ai + 02 — 1- Hence, ai = —2 
from the second equation. Substituting this into the first equation gives 02 > —2. 

Substituting (|3.2|) into (|3.18jl and requiring that at least two leading terms balance 
gives us two branches: ai = 02 = —2 and ai = —2, 02 = —1. The branches with ai = —2 
and a2>0 are excluded for they require that either no(x,t) or vo{x,t) is identically zero. 

Continuing with the two branches and solving for uo{x,t) and vo{x,t) gives 



ai = 02 = —2, 

uo{x,t) = -4gl{x,t), and 
^vo{x,t) = ±2-/6ag'^{x,t), 



ai = —2, 02 = —1, 

uoix,t) = -2gl{x,t), (3.19) 
^VQ{x,t) arbitrary. 



For the branch with ai = 02 = —2, substituting (|3.3() . 

|n(x, t) = -4.gl{x, t)g-'^{x, t) + n^(x, t)g'-'^{x, t), 
yv{x,t) = ±2\/6ag^(x, t)5~^(x, t) + Vr{x,t)g^^'^{x,t), 

into H3.18|l and equating to zero the coefficients of the lowest order terms in g{x, t) that 
are linear in and Vr gives 

-{r-A){r^ -hr-lS)agl{x,t) ±l2VQ^gl{x,t) \ (ur{x,t)\ ^ 

^A{r - A)V^gl{x,t) [r - 2)(r - 7)rgl[x,t) ) \vr{x,t) ) ' ^ 

dei{Qr) = -a{r + 2){r + l){r-Z){r-A){r-Q){r-S)gl{x,t) = 0, (3.22) 

we obtain the resonances ri = —2, r2 = — 1, T3 = 3, r4 = 4, = 6, and rg = 8. 

By convention, the resonance ri = —2 is ignored since it violates the hypothesis that 
g{x,t)~'^ is the dominant term in the expansion near g{z) = 0. Furthermore, this is not a 
principal branch since the series has only five arbitrary functions instead of the required six 
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(as the term corresponding to resonance ri = —2 does not contribute to the expansion). 
Thus, this leads to a particular solution and the general solution may still be multi-valued. 

As in the previous example, the constants of integration at level k are found by sub- 
stituting p.Sf) into H3.18() and pulling off the coefficients of g^~^{x,t). At level k = 1, 



f llagl{x,t) ±2^^Q^gl{x,t)\ fui{x,t)\ ^ f 20agl{x,t)ga:x{x,t) \ 
\±2V6^gl{x,t) -gl{x,t) J\viix,t)J \±10V6^g^^{x,t)gxxix,t)J ' ^ ■ ' 

and thus, 

ui{x,t) =4fc(x,t), vi{x,t) = ±2VQagxxix,t). (3.24) 
At level k = 2, 

^glxi^^t) - 9x{x,t){gt{x,t) +Ag3x{x,t)) 



U2{x,t) 



V2{x,t) = zb 



2>gl{x,t) 

(1 + 2a)gt{x, t)gx{x, t) + 4agx{x, t)g3x{x, t) - 3agl^{x, t) 



(3.25) 



6ag2(x,t) 



The compatibility conditions at levels A; = = 3 and k = = 4: are trivially satisfied. At 
levels k = 5 and k = 7, Uk{x,t) and Vk{x,t) are unambiguously determined (not shown). 
At resonance levels k = = 6 and k = tq = 8, the compatibility conditions require a = ^. 



(3.26) 



Likewise, for the branch with ai = —2, 02 = —1, substituting (|3.3|) . 

(u{x, t) = -2gl{x, t)g-'^{x, t) + n^(x, t)/"2(x, i), 
|t;(x,t) = vo{x,t)g^''{x,t) + Vr{x,t)g''-'^{x,t), 

into (|3.18|) gives 

-a{r + l){r - 4){r -6)g^{x,t) -3vo{x, t)gx{x,t) \ fur{x,t)\ _ , . 

r(r - l)(r - 5)gl{x,t) J \vrix,t) J ~ ^"^'^'^ 

Since 

det{Qr) = -a{r + l)r(r - l)(r - 4)(r - 5)(r - 6)g^{x,t), (3.28) 

the resonances are ri = —1, r2 = 0, = 1, r4 = 4, rs = 5, and rg = 6. 

Since r2 = is a resonance, there must be freedom at level k = r2 = 0; indeed, 
coefficient no(x,t) = —2g^{x,t) is unambiguously determined but ^0(2;,^) is arbitrary. 
Then, the constants of integration are found by substituting 

(u{x,t) = -2gl{x,t)g-^{x,t) + ui{x,t)g-'^{x,t) H ^ U(i{x,t)g^{x,t), 

\v{x,t) = VQ{x,t)g~^{x,t) + vi{x,t) H h V6{x,t)g^{x,t), 

into H3.18|) and pulling off the coefficients of g^~^{x,t) in the first equation and g^~'^[x,t) 
in the second equation. At level k = = 1, 

f a 0\/ni(x,t)\ / 2agxx{x,t) \ 
\vQ{x,t) 0j\vi{x,t)J \2vo{x,t)gxx{x,t)J ' 
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So, ui{x,t) = 2gxxix,t) and vi{x,t) is arbitrary. At level k = 2, 

f I2agl \(u2\^f 2gtgx+6aglx-Vo-8agxg3x 
\-3vogx -Qgl) \V2) \vm+^{'^i)xgl-'i{vQ)xgxx+Kvo)2xgx+VQgzx 

which unambiguously determines U2{x,t) and V2{x,t). Similarly, the coefficients in the 
Laurent series solution are unambiguously determined at level A; = 3. At resonance level 
A; = r4 = 4, the compatibility condition is trivially satisfied. At resonance levels k = = 5 
and /c = rg = 6, the compatibility conditions requires a = ^. 

Therefore, H3.18() satisfies the necessary conditions for having the Painleve property 
when fact confirmed by other analyses of complete integrability [1]. 



4 Key Algorithms 

In this section, we present the three key algorithms in greater detail. To illustrate the 
steps we consider a generalization 

iut + Uxx + {\u\'^ + /^kP)"" + o.(x, t)u + c(x, t)v = 0, 
ivt + Vxx + {\v\ + I3\u\ )v + b{x,t)v + d{x,t)u = 0, 

of the coupled nonlinear Schrodinger (NLS) equations [40] 

iut + Uxx + i\u\'^ + P\v\'^)u = 0, 
ivt + vxx + {\v\^ + f3\u\^)v = 0. 

In (|4.1j) . a{x,t), . . . ,d{x,t) are arbitrary complex functions and /3 is a real constant pa- 
rameter. Since all the functions in 1)4. 1|) are complex, we write the system as 

iut + Uxx + {uu + I3vv)u + a{x, t)u + c(x, t)v = 0, 
iut — Uxx — (uu + I3vv)u — a{x, t)u — c(x, t)v = 0, 
ivt + Vxx + {vv + f3uu)v + b{x, t)v + d{x, t)u = 0, 
ivt — Vxx — {vv + I3uu)v — b{x, t)v — d{x, t)u = 0, 

and treat ti, u, v, and v as independent complex functions. As is customary, the variables 
with overbars denote complex conjugates. 



4.1 Algorithm to determine the dominant behavior 

Determining the dominant behavior of (|2.1|) is delicate and the omission of valid dominant 
behaviors often leads to erroneous results [33]. 

Step 1 (Substitute the leading- order ansatz). To determine the values of a^, it is sufficient 
to substitute Ui{z) = ^^^(z)"', into (|2.1|) . where Xi is constant and g{z) is an analytic 
function in a neighborhood of the non-characteristic manifold defined by g{z) = 0. 

Step 2 (Collect exponents and prune non-dominant branches). The balance of exponents 
must come from different terms in ()2.1[) . For each equation = 0, collect the exponents of 
g{z). Then, remove non-dominant exponents and duplicates (that come from the same term 
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in (|2.H) ). For example, ai + 1 is non-dominant and can be removed from {ai — 1, ai + 1} 
since ai — 1 < ai + 1. 

For (|4.3() . the exponents corresponding to each equation are 



F3 
Fa 



{ai — 2, 2ai + a2, "i + "3 + "4}, 

{02 - 2,ai + 2a2,a2 + "3 + "4}, ^^^^ 
{03 — 2, 2a3 + 04, ai + 02 + 03}, 



{04 — 2, 03 + 2a4, «! + 02 + on], 
after duplicates and non-dominant exponents have been removed. 

Step 3 (Combine expressions and compute relations for ai) . For each Fi separately, equate 
all possible combinations of two elements. Then, construct relations between the ctj by 
solving for ai,a2, etc., one at a time. 
For we get 

Fi : {«! — 2 = 2ai + a2,ai — 2 = ai + a3 + 04, 

2ai + a2 = ai + as + 04} (4.5) 
^ {ai + Q2 = -2, 03 + Q4 = —2, ai + 02 = as + "4}- 

For ^2,^3 and F4 we again find that {qi -|- 02 = —2, 03 -|- 0:4 = —2, ai + 02 = + 04}. 

^tep (Combine equations and solve for exponents ai). By combining the sets of expres- 
sions in an "outer product" fashion, we generate all the possible linear equations for ai. 
Solving these linear systems, we form a set of all possible solutions for ai. 
For ()4.3() . we have three sets of linear equations 

U+a2 = -2 ^h^"^;-'' (4.6) 
[03 + a4 > -2, 

{«, + a4 = -2 ^hl^^^-l' (4.7) 

>- [03 + "4 = -2, 

and 

jai + 02 = 03 + a4, jai + Q2 = 03 + 04 > -2. (4.8) 

Although the algorithm treats u, u, v, and v as independent complex functions, we 
know that ai = a2 and a^ = 04 because u and v are the complex conjugates of u and v. 
Our package PainleveTest .m can take advantage of such additional information by using 
the option DominantBehaviorConstraints -> {alpha [1] == alpha [2] , alpha [3] == 
alpha [4] } . Using this additional information yields three cases, ai = Q2 = as = 04 > — 1 
and ai = a2 = —1, a^ = a^ > —1 and ai = 02 > —1, 03 = 04 = —1. 

Step 5 (Fix the undetermined ai). First, compute the minimum values for the unde- 
termined ai. If a minimum value cannot be determined, then the user-defined value 
DominantBehaviorMin is used. If so, the value of the free ai is counted up to a user 
defined DominantBehaviorMax. If neither of the bounds is set, the software will run the 
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test for the default values Oi = —1,-2 and —3. For maximal flexibility, with the option 
DominantBehavior one can also run the code for user-specified values of Oj. An example 
is given in Section d In any case, the selected or given dominant behaviors are checked 
for consistency with H2.1|l . 

For ()4.3() . if we take ai, . . . , 04 < 0, then we are left with only one branch 

oi = ^2 = 03 = ^4 = —1. (4-9) 

Step 6 ( Compute the first terms in the Laurent series). Using the values for aj, substitute 



ni z 



Ui,o(z)5"' (z) 



(4.10) 



into (|2.1j) and solve the resulting (typically) nonlinear equations for Uifi{z) using the 
assumption that Uifi{z) ^ 0. 
For (|0)) . we find 



ai = 02 = «3 = «4 = —1) 
uo{x,t) = -2gl{x,t){l + (3)-^UQ\x,t), 
[vo{x,t) = -2gl{x,t){l + (3r^VQHx,t), 



(4.11) 



where uo{x,t) and vo{x,t) are arbitrary functions. 

If we do not restrict ai, . . . , 04 < 0, then there are contradictions with the assumption 
Uifl{z) ^ for all but two possible dominant behaviors. 



ai = 02 
"3 > 3, 
0^4 > 3, 



-1, 



and 



ai > 3, 

02 > 3, 

03 = 0(4 



(4.12) 



-1. 



4.2 Algorithm to determine the resonances 

Step 1 ( Construct matrix Qr)- Substitute 



Ui{z) = Uifi{z)g°'' (z) + Ui^r{z)g' 



(4.13) 



into 1)2. 1|) . Then, the («,j)-th entry of the M x M matrix Qr is the coefficients of the 
linear terms in tij.r(z) of the leading terms in equation Fi = 0. 

Step 2 (Find the roots of det{Qr))- The resonances are the solutions of det{Qr) = 0. If 
any of these solutions (in a particular branch) is non-integer, then that branch of the 
algorithm terminates since it implies that some solutions of (|2.1|) have movable algebraic 
branch point. If any of the resonances are rational, then a change of variables in (|2.1jl 
may remove the algebraic branch point. Such changes are not carried out automatically. 
For branch (14.111). 



det(Q,) = (r - 4)(r - 3)^ {r\l + p) - 3r(l + /3) - 4(1 - /3)} r\r + 1) 



X {1 + pf4{x,t)v^oix,t)gUx,t). 



(4.14) 
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Since the roots of (|4.14j) for r depend on the constant parameter (3, we must choose values 
of P so that all the solutions are integers before proceeding. For (3 = 1, the resonances are 
n = -1, T2 = = ?'4 = 0, rs = rg = ry = 3, = 4. 

While taking /? = also yields all integer resonances, it violates the assumption that 
all the parameters in H2.1() are nonzero. Allowing the parameters in ()2.1|) to be zero could 
cause a false balance in Algorithm l4.11 Thus, (|2.1j) with /? = should be treated separately. 
In this example however, setting /? = does not affect the dominant behavior and the 
resonances are ri = r2 = — 1, = r4 = 0, rs = rg = 3, and ry = rg = 4. 

Although taking (3 = 25/7 leads to rational resonances at r4 = rs = 3/2, they are not 
easily resolved by a change of variables in (|4.3j) . The branches with dominant behavior, 
ai = 02 = —1, as > 04 > 3, have resonances n = — 03 — 1, r2 = — Q3+2, r-^ = —a^ — l, = 
—04 + 2,r5 = —l,rQ = 0,rj = 3 and rg = 4. Since ri,r2 < —1 when 03 = 04 = 3, 
ri,r2,r3 < —1 when 03 > 04 = 3, and ri,...,r4 < —1 when 03,04 > 3, these are not 
principal branches and should be investigated using the perturbative Painleve approach [9]. 

4.3 Algorithm to determine the constants of integration and check com- 
patibihty conditions 

Step 1 ( Generate the system for the coefficients of the Laurent series at level k ). Substitute 



into (|2.H) and multiply Fi by g~'^'-{z), where 7^ is the lowest exponent of g{z) in Fi. The 
equations for determining the coefficients of the Laurent series at level k then arise by 
equating to zero the coefficients of g^{z). These equations, at level k, are linear in Ui^k{z) 
and depend only on Uij{z) and g{z) (and their derivatives) for 1 < i < M and < j < k. 
Thus, the system can be written as 



where u^. = (ui^fc(z), . . . , MM,fc(z))'^- 

Step 2 (Solve the linear system for the coefficients of the Laurent series). If the rank of 
Qk equals the rank of the augmented matrix {Qk\Gk), solve (|4.16|) for the coefficients of 
the Laurent series. If A; = rj, check that rankQfc = M — Sj, where sj is the algebraic 
multiplicity of the resonance rj in det{Qr) = 0. 

If rankQfc 7^ rank((5fc|Gfc), Gauss reduce the augmented matrix {Qk\Gk) to determine 
the compatibility condition. If all the compatibility conditions can be resolved by restrict- 
ing the coefficients parameterizing 1)2. 1() . then ()2.1|) has the Painleve property for those 
specific values. If any of the compatibility conditions cannot be resolved by restricting 
the coefficients parameterizing ()2.ip . then the Laurent series solution for this branch has 
a movable logarithmic branch point and the algorithm terminates. 




(4.15) 



fc=0 



QfcUfe = Gfc(uo, ui, . . . , Uk-i,g, 



(4.16) 
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For (|4.3|) with f3 = 1, the principal branch 



ai = Q2 = as = Q!4 = —1, 

uo{x,t) = -u^^{x,t){vo{x,t)vo{x,t) - 2gl{x,t)}, 

uo{x, t), vo{x, t),vo{x, t) arbitrary, 

n = -1,-^2 = ?"3 = ^'4 = 0, rs = re = ry = 3, rs = 4, 



(4.17) 



has three compatibility conditions at level k = = = r-j = 3. These compatibility 
conditions require that ax{x,t) = ax{x,t) = bx{x,t) = bx{x,t) and Cx{x,t) = dx{x,t) = 0. 
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At level k = rg = 4, the compatibility condition requires d{t) = c{t) and 

{{a - af + 2i(o^ - a^)h'{t)}{2 + voVq) - {{b - If + 2i(h^ - h)h'{t)}voVo 

+ 2i{a -a + b- b){vo{vo)t + ivo)tVo) + i{at - at - h + bt)voVo (4.18) 

+ 2i{at - at) - (ar^x + a^x - bxx - bxx)voVo - 2axx + Qaxx = 0, 

where we have taken g{x,t) = x — h{t). Careful inspection of H4.18() reveals that a{x,t) = 
b{x,t). Setting a{x,t) = b{x,t) = r{x,t) + is{x,t), where r{x,t) and s{x,t) are arbitrary 
real functions, H4.18() becomes 

2s'^{x, t) + st{x, t) + 2h'{t)sx[x, t) - Txxix, t) + 2isxx{x, t) = 0. (4.19) 

Since h'{t) is arbitrary, it follows that Sx{x,t) = 0. Thus, rxx{x,t) = 2s^{t) + s'{t) and 
upon integration 

r{x,t) = ]^{2s^{t)+s'{t)}x'^ +ri{t)x + r2{t), (4.20) 

where ri[t) and r2{t) are arbitrary functions. 

Therefore, the generalized coupled NLS equations, 

iut + Uxx + (|uP + \v\^)u + {{s^{t) + \s{t)]x^ + ri{t)x + r2{t) + is{t)] u + c{t)v = 0, 

ivt + Vxx + {\u\^ + \v\^)v + {{s^{t) + \s'{t)]x'^ + ri{t)x + r2{t) + is{t)} v + c{t)u = 0, 

passes the Painleve test, where ri(t),r2(t), and s{t) are arbitrary real functions and c{t) 
is an arbitrary complex function. 

When /3 = 0, the two compatibility conditions at level k = = = 2> require that 
c(x, t) = (i(x, t) = c(x, t) = d{x^ t) = 0. Similarly, the compatibility conditions at level 
= ry = rg = 4, require that 

a{x, t) = {s^{t) - y{t)]x'^ + ri{t)x + r2{t) + is{t), (4.21) 
where ri{t),r2{t) and s{t) are arbitrary real functions. Therefore, 

iut + Uxx + \u\'^u + {{s^{t) - y{t)]x^ + ri{t)x + r2{t) + is{t)} n = 0, (4.22) 
passes the Painleve test, a fact confirmed in [1]. 



5 Additional Examples 

5.1 A peculiar ODE 

Consider the ODE [34] 

u^u'" - ^{u'f = 0. (5.1) 

Substituting (|XT|) into gives a{a + 2){2a-l)x^g{zf'^°'^^'> = 0. So, both the terms in 
1)5. 1|) have the same leading exponent, 3(a — 1). Using the procedure in Section in Step 
5 the software automatically runs the test for the default values a = —3, —2, and —1. The 
choices q = — 1 and —3 are incompatible with the assumption uq 7^ 0. The leading term 
vanishes for a = —2 and uq is arbitrary. Substituting u{z) = uog~'^{z) + ui^r 9^~'^{z), we 
find that ri = — l,r2 = 0, and = 10. Thus, the Laurent series solution of (|5.1|) is 

u{z) = uo{z - zo)"^ + uioiz - zof H , (5.2) 

where zq,uq and uiq are arbitrary constants. Hence, ()5.1|) passes the Painleve test. 
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5.2 The sine-Gordon equation 

Consider the sine-Gordon equation [1], 

uu + Uxx = smu. (5.3) 
Using the transformation v{x,t) = e*"(^'*), we obtain a polynomial differential equation 

vvtt + vvxx-v^ -vl = ]^v{v'^ -I). (5.4) 

The dominant behavior of (|5.4|) is v{x,t) ~ 'i{g'^{x,t) + g^{x,t))g'^'^{x,t), with resonances 
ri = — 1 and r2 = 2. The Laurent series solution of 1)5. 4|1 is 

V = A{gl + gl)g~^ - iig^x + gtt)g'^ + ^2 + • • • , (5.5) 

where g{x,t) and V2{x,t) are arbitrary functions. The sine-Gordon equation passes the 
Painleve test and is indeed completely integrable [1]. 

5.3 The cyUndrical Korteweg-de Vries equation 

Consider the generalized KdV equation, 

ut + 6uux + + a{t)u = 0, (5-6) 

where a{t) is an arbitrary function parameterizing the equation. The dominant behavior 
of 1)5. 6() is u{x,t) ~ —2g'^{x,t)g~'^{x,t), with resonances ri = — l,r2 = 4 and = 6. At 
level A; = r3 = 6, we obtain the compatibility condition 

^''<"'+°'">^0. (5.7) 

Qgx[x,t) 

So, 1)5. 6|) passes the Painleve test if a{t) = ^. In this case, (|5.6() reduces to the cylindrical 
KdV, which is completely integrable as confirmed by other analyses [1]. 

5.4 A fifth-order generahzed Korteweg-de Vries equation 

Consider the generalized fifth-order KdV equation, 

Ut + aUxUxx + buU3x + CU^Ux + U5x = 0, (5.8) 
with constant parameters a, b, and c. The dominant behavior of (|5.8j) is 

u{x,t) ~ {(a + 2b) ± + ^ab + - 40c} g-^{x, t). (5.9) 

The resonances are the roots of 



det(Q^) = -c(r - 6)(r + 1) (^3\/(a + 26)2 - 40c(2a - 6(r - 4)) 

- 6(a + 26)2 + 240c + (36(a + 26) - 86c)r + IScr^ - cr^) g^. (5.10) 
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Determining what values of o, b, and c that lead to integer roots of (|5.1())) is difficult by 
hand or with a computer. An investigation of the scaling properties of ()5.8() reveals that 
only the ratios a/b and are important. Let us consider the well-known special cases. 

If we take a = b and 5c = 6^, then (|5.8() passes the Painleve test with resonances 
n = -2,r2 = -l,r3 = 5,r4 = 6, rs = 12 and n = -l,r2 = 2,r3 = 3, r4 = 6,r5 = 10. 
Taking 5 = 5, equation (|5.8() becomes the completely integrable equation 

ut + 5uxUxx + 5nu3^ + 5n^na; + ttsx = 0, (5-11) 

due to Sawada and Kotera [36] and Caudrey et al. [5]. 
If we take a = 2b and 10c = 36^, then (gSl) 

passes the Painleve test with resonances 
ri = — 3,r2 = —1,^3 = 6,r4 = 8,r5 = 10 and ri = — l,r2 = 2,r3 = 5,r4 = 6,r5 = 8. For 
6 = 10, equation (|5.8() is a member of the completely integrable KdV hierarchy 

Ut + lOuU^x + 20UxUxx + SOn^Ua; + tiSa; = 0, (5-12) 

due to Lax [26] . 

If we take 2a = 55 and 5c = b^, then (|5.8() passes the Painleve test with resonances 
ri = -7, r2 = -l,r3 = 6,r4 = 10, rs = 12 and ri = -l,r2 = 3, r3 = 5,r4 = 6,r5 = 7. 
When b = 10, equation (|5.8() is the Kaup-Kupershmidt equation [12,17], 

Ut + 10uU3a. + 2f)u^Ux + 25na.lia:a: + U^x = 0, (5.13) 

which is also known to be completely integrable. 

While there are many other values for a, 6, and c, for which 1)5. 10() only has integer roots, 
but compatibility conditions prevent (|5.8|) from having the Painleve property. For instance, 
when a = 2b and 5c = 26^, the resonances are ri = — l,r2 = 0, r3 = 6, r4 = 7, rs = 8. 
At level A; = r2 = 0, we are forced to take UQ{x,t) = —30g'^{x,t)/b, so the Laurent series 
solution is not the general solution and (|5.8|) fails the Painleve test. Similarly, when 
7a = 19b and 49c = 96^, we have resonances ri = — 1, r2 = 3 and = = = 6, so the 
Laurent series solution is not the general solution and, again (|5.8() fails the Painleve test. 

6 Brief Review of Symbolic Algorithms and Software 

There is a variety of methods for testing nonlinear ODEs and PDEs for the Painleve prop- 
erty. While the WTC algorithm discussed in this paper is the most common method used 
in Painleve analysis, it is not appropriate in all cases. For instance, there are numerous 
completely integrable differential equations which have algebraic branching in their series 
solutions; a property that is allowed by the so-called "weak" Painleve test (see [13,32,33]). 
A more thorough approach for testing differential equations with branch points is the poly- 
Painleve test (see [22,23]). The perturbative Painleve test [9] was developed to check the 
compatibility conditions of negative resonances other than r = — 1. 

For testing ODEs, there are several implementations: ODEPAINLEVE developed by Rand 
and Winternitz [34] in Macsyma is restricted to scalar differential equations; PTEST.RED 
by Renner in Reduce [35]; and, a Reduce package by Scheen [37] which implements both 
the traditional and the perturbative Painleve tests. For testing PDEs, there are a few 
implementations. The package PAINMATH . M by Hereman et al. [16] is unable to find all the 
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dominant behaviors in systems with undetermined and is limited to two independent 
variables. 

Only the Maple package PDEPtest by Xu and Li [45-47] is comparable to our package 
PainleveTest .m [4]. The package PDEPtest was written after our package and allows 
the testing of systems of PDEs (but not ODEs) parameterized by arbitrary functions us- 
ing either the traditional WTC algorithm or the simplification proposed by Kruskal (see 
Sectional). While PDEPtest can find all the dominant behaviors in some systems with 
undetermined Oi (such as the Hirota-Satsuma system), it fails to find the dominant be- 
haviors for systems in which more than one is undetermined (such as the NLS equation, 
iut+Uzz + '^u\u\'^ = 0, which is completely integrable [1]). Furthermore, PDEPtest requires 
that all the ai are negative, a weakness of the implementation, since it is standard to allow 
some positive exponents (see equation (2.4) in [33] with leading exponents —1 and 1). 



7 Using the Software Package PainleveTest .m 

The package PainleveTest .m has been tested on both PCs and UNIX work stations with 
Mathematica versions 3.0, 4.0, 4.1, 5.0, 5.1, and 6.0 using a test set of over 50 PDEs and 
two dozen ODEs. The Backus-Naur form of the function is 

{Main Function) — > PaLinleveTest[{Equations), {Functions), 

(Variables) , {Options)] 

{Options) — > Verbose {Boolean) \ 

KruskalSimplif ication ^ {Variable) \ 
DominantBehaviorMin {Negative Integer) \ 
DominantBehaviorMax — > {Integer) \ 
DominantBehavior — > {List of Rules) \ 
DominantBehaviorConstraints ^ {List of Constraints) \ 
DominantBehaviorVerbose {Range) \ 
ResonancesVerbose {Range) \ 
ConstantsOf IntegrationVerbose — > {Range) 

{Bool) True ] False 

{Range) ] 1 ] 2 ] 3 

{List of Rules) {{alpha[l] ^ (Integer ), alpha[2] -^{Integer), ...}, ...} 

{List of Constraints) — > {alpha[l] == alpha[2], ...} 



The output of the function is 

I {{Dominant behavior}, {Resonances}, 

{{Laurent series coefficients}, {Compatibility conditions}}, . . . | 

If using a PC, place the package PainleveTest .m in a directory, say myDirectory on drive 
C. Start a Mathematica notebook session and execute the commands: 



In[l] = SetDirectory ["c : WmyDirectory"] ; (* Specify the directory *) 
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In [2] = Get ["PainleveTest .m"] (* Read in the package *) 

In [3] = PainleveTest [ (* Test the KdV equation *) 

{D[u[x,t] ,t]+6*u[x,t]*D[u[x,t] ,x]+D[u[x,t] ,{x,3}] == 0}, 
u[x, t] , {x,t}, KruskalSimplif ication -> x] 

Out [3] = 

{{{ai ^ -2},{r ^ -l,r ^ 4,r ^ 6}, 

h'(t) 

{{Ulfl -2,Mi,i 0,Ui,2 —^,ui^3 0, 



ni,4 ^ Ci(t),ni,5 - ^,^^1,6 - C2{t)},{}}}} 

The option KruskalSimplif ication -> x allows one to use g{x,t) = x — h{t) in the 
calculation of the constants of integration and in checking the compatibility conditions. 

In [4] = PainleveTest [ (* Eq. (2.4) in Ramani et al. [32] *) 
{D[x[z], z] == x[z]*(a - x[z] - y[z]), 

D[y[z] , z] == y[z]*(x[z] - 1)>, 
{x[z], y[z]}, {z}, DominantBehaviorMax -> 1 ] 

Out [4] = 

{{{ai ^ -1,02 ^ -l},{r ^ -l,r ^2},{{ni,o ^ -1, . . . }, {a + 1 = 0}}}, 
{{qi ^ -1,02 ^ l},{r ^ -1,^ 0},{{ni,o ^ l,U2,o ^ Ci},{}}}} 

In this example, if the DominantBehaviorMax option was not used, we would wrongly 
conclude that the system only passes the Painleve test when a = —1. However, by allowing 
positive ai, we find the second branch oi = — 1 and a2 = 1, for which the system passes 
the Painleve test without restricting the value of the parameter a. Alternatively, executing 

In [5] = PainleveTest [{ D[x[z], z] == x[z]*(a - x[z] - y[z]), 
D[y[z], z] == y[z]*(x[z] - 1)}, {x[z], y[z]>, {z>, 
DominantBehavior -> {{alpha [1] -> -1, alpha [2] -> 1}} ] 

would only test the branch with oi = — 1 and 02 = 1- For an example of the option 
DominantBehaviorConstraints, see Step 4 of Algorithm 14.11 

The option Verbose -> True gives a brief trace of the calculations in each of the three 
steps of the algorithm. The options DominantBehaviorVerbose, ResonancesVerbose, 
and ConstantsOf IntegrationVerbose allow for a more detailed trace of the calcula- 
tion. For instance, DominantBehaviorVerbose -> 1 would show the result of substi- 
tuting the ansatz, the exponents before and after removing non-dominant powers, etc. 
While DominantBehaviorVerbose -> 3 shows the result of nearly every line of code in 
the package, allowing the user to check the results in the trickiest cases. 
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8 Discussion and Conclusions 

Our software package PainleveTest .m is applicable to polynomial systems of nonlinear 
ODEs and PDEs. While the Painleve test does not guarantee complete integrability, it 
helps in identifying candidate differential equations for complete integrability in a straight- 
forward manner. For differential equations with parameters (including arbitrary functions 
of the independent variables), our software allows the user to determine the conditions 
under which the differential equations may possess the Painleve property. Therefore, by 
finding the compatibility conditions, classes of parameterized differential equations can be 
analyzed and candidates for complete integrability can be identified. 

The difficulty in completely automating the Painleve test lies in determining the domi- 
nant behaviors of the Laurent series solutions; specifically, determining all the valid domi- 
nant behaviors when one or more of the Oi are undetermined. While there are other imple- 
mentations for the Painleve test, ours is currently the only implementation in Mathematica 
which allows the testing of polynomial systems of nonlinear PDEs with no limitations on 
the number of differential equations or the number of independent variables (except where 
limited by memory). 
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